Mapping Itinerant Electrons around Kondo Impurities 
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We investigate single Fe and Co atoms buried below a Cu(lOO) surface using low temperature 
scanning tunneling spectroscopy. By mapping the local density of states of the itinerant electrons 
at the surface, the Kondo resonance near the Fermi energy is analyzed. Probing bulk impurities 
in this well-defined scattering geometry allows separating the physics of the Kondo system and the 
measuring process. The line shape of the Kondo signature shows an oscillatory behavior as a function 
of depth of the impurity as well as a function of lateral distance. The oscillation period along the 
different directions reveals that the spectral function of the itinerant electrons is anisotropic. 
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For a long time it has been well known that a local- 
ized spin degree of freedom — a magnetic impurity — in a 
nonmagnetic host metal significantly alters the scattering 
behavior of the conduction band electrons of the host as 
compared to nonmagnetic impurities. This results in a 
variety of thermodynamic anomalies, which are summa- 
rized by the term Kondo effect [1 . The most prominent 
macroscopic hallmark is the resistance minimum at low 
temperatures observed for metals with magnetic impu- 
rities. From a microscopic point of view the impurity 
interacts with the surrounding electron gas — the itiner- 
ant electrons. Below a characteristic temperature, the 
Kondo temperature Tk, a narrow many-body resonance 
named Kondo or Abrikosov-Suhl resonance builds up in 
the spectral function of the impurity at the chemical po- 
tential and the impurity spin is effectively screened. Elec- 
tric transport is dominated by electrons near the Fermi 
energy for low temperature. Hence the Kondo resonance 
leads to a strong scattering of electrons at the impurity 
and an increase of resistivity with decreasing tempera- 
ture. 

Microscopic properties of single impurity Kondo sys- 
tems on the atomic scale regained interest by recent scan- 
ning tunneling microscopy (STM) experiments [2-7 on 
single magnetic ad-atoms and molecules deposited on no- 
ble metal surfaces (for a review see [8 ). In these works 
the Kondo effect manifests itself as a sharp signature 
in the STM differential conductance around zero bias. 
Studies which investigate the dependence of the Kondo 
signature on the lateral distance [2, 5^ showed that the 
signal rapidly vanishes and the line shape nearly remains 
constant. This is in contrast to theory [9HTT] which pre- 
dicts a long range visibility and oscillatory behavior of 
the spectral function of the itinerant electrons - the local 
density of states (LDOS). Models for tunneling through 
ad- atom systems have to treat the tip, its coupling to the 
localized impurity state, its coupling to the itinerant bulk 
electrons as well as surface states, and all interferences 



between alternative transport paths [TOl [TT] . 

During recent years STM was refined to investigate 
nano structures p!2UT6i or even single impurity atoms 
pm2Q] which are buried below a metal surface. The 
interference pattern on the surface caused by the sub 
surface scattering centers are determined by the band 
structure of the underlying host material [17^. This in- 
duces strong directional Friedel oscillations with wave 
vectors and amplitudes depending on the crystal direc- 
tion [20\. Recently, it was shown that for magnetic im- 
purities buried below the surface a long range signature 
of the Kondo effect is observable [21]. 

In this letter we present a detailed analysis of differ- 
ential conductance as a function of lateral distance for 
Fe and Co atoms situated at different depths below a 
Cu(lOO) surface. We show that the distance dependence 
of the line shape of the Kondo signature is directly con- 
nected to the band structure of the host crystal. In con- 
trast to previous works using ad-atom systems which had 
studied only local properties, this work focuses on sin- 
gle Kondo impurities in the bulk on length scales, which 
hitherto were not accessible. The long range Kondo sig- 
nature allows us to explore the full microscopic properties 
around magnetic impurities. Additionally, probing bulk 
impurities with STM allows separating the physics of the 
correlated many-body system and the measuring process. 

The experiments were performed using a home-built 
low-temperature scanning tunneling microscope (STM) 
operating at 6K at pressures below 5 • 10~^^ mbar. The 
Cu(lOO) single crystal substrate is cleaned by several cy- 
cles of argon bombardment and electron beam heating. 
A copper alloy with a small amount (~ 0.02%) of im- 
purities is prepared through simultaneous deposition of 
copper and the impurity material from two electron beam 
evaporators. 

The embedded atoms are identified by their short pe- 
riod standing-wave pattern which has four fold symme- 
try [Fig. [Ik]. The modulation in the topographic images 
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FIG. 1: (a) Scheme of the experiment: the green circle marks 
the impurity position below the surface and constant current 
topographies (2nm x 2nm, V = lOmV, I = 2nA) of single 
Co atoms buried below a Cu(lOO) surface, (b) Corresponding 
cross sections (shifted vertically for clarity) of each layer along 
the X = [010] direction, marked by the arrows in (a). 



has a wavelength of 5.5 ± 0.4 A and the amplitude is of 
the order of a few pm [see Fig.[l]3]. Since the surface 
is atomically flat and the applied bias voltage is small, 
the modulation amplitude in the topographic images can 
be directly connected to a change in the local density 
of states LDOS. Assuming a work function of 4.3eV, a 
change of 1pm in tip-sample distance results in a change 
in the LDOS of 2%. Consequently, the LDOS modu- 
lations that cause the observed modulation amplitudes 
are in the range of 6% (7ML) - 13% (3ML). The depth 
(in expressed in monolayer (ML) of each impurity be- 
low the surface is determined as follows: the observed 
standing wave patterns are ordered by the lateral size. 
The size of the standing wave pattern increases mono- 
tonically with depth of the impurity. As a verification 
of the obtained impurity position we use atomically re- 
solved topographies. Since the Cu(lOO) has an fee crystal 
structure and impurity atoms are located on lattice sites, 
the impurity contrast follows a certain pattern with re- 
spect to the host lattice. We label the surface layer as 1st 
monolayer. The standing wave pattern of an odd layer 
impurity is centered on a corrugation maximum of the 
surface. If the impurity is positioned in an even layer the 
center is located between the corrugation maxima. 

To investigate the embedded impurities and their influ- 
ence on the LDOS at the surface we use scanning tunnel- 
ing spectroscopy (STS) data by recording an I{V) curve 
at every scanning point with interrupted feedback loop. 
Further data processing includes averaging and numeri- 



cal differentiation. This provides a complete map of the 
differential conductance dl /dV{x,y^V) as a function of 
lateral tip position (x, y) on the surface and applied bias 
voltage V. In order to perform a quantitative analysis of 
the experimental data, drift effects which occur during 
acquisition of spectroscopic data are removed by com- 
paring the simultaneously measured topography with a 
reference topography obtained in the normal constant 
current mode. 
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FIG. 2: (a) Spectroscopy of a 3ML Fe impurity, the cor- 
responding topography (2.5nm x 2.0nm, V = -312mV, I = 
1.8nA) is shown in the inset. A pair of dl/dV spectra taken 
with the STM tip held over the center (purple curve) of the in- 
terference pattern, evaluated at the purple cross and over the 
nearby bare Cu surface (red curve) evaluated above the undis- 
turbed region of the image (red rectangle), (b) Difference 
of the differential conductance AdI/dV = dl/dV — dIo/dV 
(green curve) above the impurity (dl/dV) and above the free 
surface (dlo/dV). The solid black curve shows a fit to the 
AdI/dV data described in the text. 

Typical differential conductance dl/dV spectra mea- 
sured directly at the center of the interference pattern 
reveal a strong signature around zero bias voltage which 
is attributed to the Kondo effect and cannot be seen in 
spectra of the bare Cu surface [see Fig.[2|i]. To extract 
the LDOS change and to remove artifacts which originate 
from the bare Cu(lOO) surface and the tip, we normalize 
the obtained differential conductance around an impurity 
by subtracting the spectra of the free surface d/o / d^ far 
away from the impurity. The difference is proportional 
to the change in the LDOS at the lateral tip position [22] . 
The normalized AdI/dV spectrum [Fig.[2]3] only shows 
a dip around zero bias, all other signatures coming from 
the free surface and tip are removed. 

Although the Kondo effect has been subject of a wealth 
of experimental and theoretical investigations, micro- 
scopic real space descriptions are rare. The most relevant 
property is the spectral function of the localized orbital, 
which has been directly measured by high-resolution 
photo emission electron spectroscopy and inverse photo 
emission [23]. In an STM experiment the current is deter- 
mined by tunneling into the conduction electron density 
of states. Therefore, the STM is probing the spectral 
function of the itinerant electrons and not the localized 
ones. A first theoretical description of the Kondo effect 
in the spectral function of the itinerant electrons at dif- 
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ferent distances from the impurity was given by [9HTT], 
more recent results can be found in [24l |25]. They all 
showed that the Kondo resonance leads to a scattering 
resonance in the LDOS resulting in an enhanced scatter- 
ing amplitude and an energy-dependent phase shift at the 
Fermi energy [21]. Here we show, that the Kondo signa- 
ture shows an oscillatory behavior as a function of depth 
of the impurity as well as function of lateral distance. 
Furthermore we analyze the influence of the electronic 
band structure of the host material on the line shape 
of the Kondo signature. As described in [9] the LDOS 
at different positions depends on the propagation of the 
conduction band electrons between tip and impurity (a 
property of the host metal) and the scattering behav- 
ior of the impurity, which is strongly influenced by the 
Kondo resonance. We replace the Lorentzian shape of 
the Kondo resonance by a phenomenological form found 
by Frota et al. [26l|27]. Their formula gives a better de- 
scription of STM experiments [2TJ^ and is justified by 
theoretical many-body numerical renormalization group 
(NRG) calculations [29]. The final expression for Frota's 
phenomenological line shape has the following functional 
form: 
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Here ek is the position of the Kondo resonance. The reso- 
nance width F is proportional to the Kondo temperature 
Tk, whereas the exact prefactor is still an open issue. 
We use F = 1.43 /cB^k, which results from NRG calcu- 
lations [30]. A linear voltage slope is added in Eq. ([T]) 
to account for additional and approximately energy inde- 
pendent background scattering processes. Measurements 
taken on different impurity depth and spatial positions 
yield a resonance width F(Fe) = 4 ± 2meV and ek(Fe) = 
± 3meV for Fe impurities and F(Co) = 81 ± IQmeV 
and ek(Co) = —65 ± 25meV for Co atoms. Our obser- 
vation that for Co the Kondo resonance is situated well 
below the Fermi energy is in good agreement with recent 
calculations |31) . 

In the following we will show that for a magnetic atom 
the line shape of the dl/dV spectra strongly depends on 
the depth of the impurity as well on the lateral distance 
[32] . The line shape is determined by (j) in Eq. ([T]). The 
line-shape parameter (LSP) results for even multiples of 
TT {(j) = 0,27r, ...) in a dip and for odd multiples of tt 
{(j) = TT, 37r, ...) in a peak. 

First, we analyze the dependence on the lateral posi- 
tion X (the [010] direction, see sketch in Fig.[l^) while 
keeping the depth of the impurity fixed. To illus- 
trate this, the topographies of a 7th Co and 4th layer Fe 
impurity are shown, respectively in Fig. [3^ and Fig.|3]3. 
A series of Ad//dV spectra taken with the STM tip 
held at various lateral spacing from the center of the in- 
terference pattern is shown in Fig.|3j3 and Fig.lsji. At 
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FIG. 3: Line-shape variation as a function of the lateral tip 
position, (a) STM constant current topography of a 7th layer 
Co impurity (3nm x 3nm, V = -300mV, I = 2nA) and (b) 
a 4th monolayer Fe impurity (2nm x 2nm, V = -150mV, I 
= 2nA), in which the tip position for the different spectra 
are depicted. Single spectra, displaced vertically for clarity, 
for different lateral distances, measured from the center along 
the [010] direction (c) for the 7ML Co impurity and (d) for 
the 4ML Fe impurity. The solid black curve shows fits to 
the AdI/dV described in the text. Spatial dependence of the 
line-shape parameter (LSP) for 7ML Co (e) and for 4ML Fe 
(f). Error bars are obtained by the confidence interval of the 
non linear fit. The dashed lines show a model presented in 
the text. 



the center (x, ^ = 0) directly above a 4th layer impu- 
rity [Fig.|3]i] a peak is observed, which corresponds to 
a value of ^ = 0.837r. Moving away from the center the 
line shape becomes more asymmetric, turns into a nega- 
tive peak and back into a symmetric peak at a distance 
of X = 6.4A. The line shape parameter increases con- 
tinuously with the lateral position [Fig.[3]F]. The same 
behavior is observed for a Co atom [see Fig. [3^] . 

The line shape parameter should depend linearly [9] 
on the distance r to the impurity with a slope given by 
2 times the Fermi wave vector kp, i.e., (r) = 2 /cp r. 
The distance r = ^/d'^-\- x'^ in our experiment is deter- 
mined directly from the lateral distance x and the depth 
dn of the impurity below the surface. We use the follow- 
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ing relation ^ = 2 /cp ^yd'^-\- x'^ + 0o to fit the line-shape 
parameter. The offset accounts for the line shape 
directly at the impurity position. For a 4th layer Fe im- 
purity ((^4 = 5. 4 A) the best fit is obtained with param- 
eters (j)o = -0.2 ± 0.47r and = 0.98 ± O.OSA-^ For 
a 7ML Co atom {d^ = 10. 8A) the best parameters are 
^0 = 0.1 ± O.Stt and = 0.85 ± O.ISA-^. 

All theoretical models so far describing the line-shape 
dependence as a function of distance treat a homogeneous 
and isotropic electron gas model with a direction inde- 
pendent kp. However, band structure calculations for Cu 
yield values of k^ = LSA'^ in [110] and k^ = 1.45A-^ 
in [100] direction The difference between our ex- 

perimental and the theoretical value is not only observed 
for the two data sets shown here, but for all our exper- 
imental findings [34]. The discrepancy can be basically 
explained by two additional ingredients: the anisotropic 
band structure of copper, i.e., the focusing effect and the 
imaging process of the STM. The noncollinear direction 
of k-p and the group velocity results in an effective re- 
duction of the observed oscillation wave vector [20]. In 
addition, not all wave functions contribute equally to the 
tunnel current. Especially states with higher k\\ compo- 
nents decay faster into the vacuum. This implies that 
for larger lateral distances of the tip states with different 
|/cf| contribute to the tunnel current. 
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FIG. 4: Experimentally determined line-shape parameter 
(LSP) for Co atoms as function of the impurity position dn 
below the surface (the corresponding topographies are shown 
in Fig.[T|). Measured AdI/dV (green) and calculated spectra 
(black) with a tip position fixed at the center of the inter- 
ference pattern. Inset: The obtained LSP (purple dots) de- 
creases linearly. Error bars are obtained by the confidence 
interval of the non linear fit. 



A way to analyze the line-shape dependence in a sim- 
pler geometry is to look at spectra with different impurity 
depth dn while the lateral tip position is kept constant at 
the center {x^y = 0) of the interference patterns [Fig. [4]. 
In this geometry the contribution of tunneling paths is 
the same for all observed spectra. The spectrum for a 
3ML Co atom reveals an asymmetric dip. For deeper 
impurity positions the line shape becomes more asym- 
metric returning into an asymmetric dip for the 7ML, 
seen before for the 3ML, indicating an oscillation length 
of nearly 4ML. The line-shape parameter [inset of Fig. [4] 
apparently decreases linearly with increasing depth. A 
linear fit reveals a slope of -0.45 ± 0.027r/ML. This 
negative slope is due to aliasing between the period- 
icity of the lattice and the wave length of the oscil- 
lations. The line shape parameter is 27r periodic and 
hence the observations are also equivalent to a slope of 
1.55 ± 0.027r/ML, which results in an Fermi wave vector 
of /cf = 1.35 ± 0.05A~^. This is in good agreement with 
the expected value of kp = 1.45A~^ for copper in [100] 
direction [33] . 

Finally, we would like to discuss our results in compar- 
ison to previous results on ad-atom systems. In the de- 
scription of ad-atom spectra the line shape of the Kondo 
resonance is discussed in terms of interfering channels 
considering tunneling into the localized orbital and the 
conduction band states, respectively. Thus, the line 
shape is a feature of the measurement process and is as- 
sociated with the relative strength of the hybridization 
between the tip, sample and impurity. In contrast, the 
investigation of bulk impurities allows a clear spatial sep- 
aration between the Kondo impurity and STM process. 
Within our interpretation, the STM is a noninvasive tool 
that exclusively probes the LDOS of the conduction band 
electrons — the itinerant electrons. 

In conclusion, our results show that STM is able to 
characterize bulk impurities on relatively long distances 
from the impurity, which opens a new way to Kondo 
physics in real space. Using STS we have demonstrated 
that single magnetic Fe and Co impurities embedded be- 
low the Cu(lOO) surface show a complex spatially depen- 
dent Kondo signature. The spatial distribution of the 
Kondo signature is — in addition to the Kondo physics — 
strongly infiuenced by the anisotropy of the band struc- 
ture of the host crystal. This work was supported by 
the Deutsche Forschungsgemeinschaft through SFB 602 
Project A3. 
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